UCSD/PTH/98-36 

The glueball spectrum from an anisotropic lattice study HLRZ1998/62 

Colin J. Morningstar 
Dept. of Physics, University of California at San Diego, La Jolla, California 92093-0319 

Mike Peardon 

John von Neumann-Institut fur Computing and DESY, Forschungszentrum Jiilich, D-52425 Jiilich, Germany 

(January 7, 1999) 

The spectrum of glueballs below 4 GeV in the SU(3) pure-gauge theory is investigated using 
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states from two-glueball and torelon-pair states. Our determination of the spectrum significantly 
improves upon previous Wilson action calculations. 
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I. INTRODUCTION 

The gluon self-coupling in quantum chromodynamics 
(QCD) suggests the existence of glueballs, bound states 
of mainly gluons. Incontrovertible experimental evidence 
for their existence remains elusive, however. A primary 
reason for this is the difficulty in extracting the proper- 
ties of glueballs from the QCD lagrangian. Investigating 
glueball physics requires an intimate knowledge of the 
confining QCD vacuum, and such knowledge cannot be 
obtained using standard perturbative techniques. Nu- 
merical simulations of the theory on a space-time lattice 
currently provide the most reliable means of studying 
glueballs. However, correlation functions of gluonic exi- 
tations are notoriously difficult quantities to measure in 
Monte-Carlo simulations, requiring large-scale computer 
resources when applying standard stochastic techniques. 
Recently, the use of spatially coarse, temporally fine lat- 
tices and improved actions was demonstrated to dramat- 
ically increase the efficiency of glueball simulations [Q . 

The objective in this paper is to apply the techniques 
of Ref. [jlj to substantially improve our knowledge of the 
glueball spectrum in the pure SU(3) gauge theory. De- 
tailed information on this spectrum is important for val- 
idating models of confined gluons and may help focus 
experimental searches for candidate glueball resonances. 
We also view this calculation as a necessary first step 
before attempting to include the effects of quarks. Note 
that unlike the quenched approximation for mesons and 
baryons, the pure glue theory is a physical quantum field 
theory with a unitary S-matrix. First, we perform six 
simulations for spatial lattice spacings ranging from 0.1 
to 0.4 fm to determine the energies of the lowest-lying 
stationary states in all of the symmetry channels allowed 
on a cubic lattice. In many channels, we also determine 



the energies of the first-excited states. Our goal is then 
to extract the masses of as many low-lying glueballs as 
possible from the 141 measurements which were made. 
Since the spectrum of glue defined in a box with peri- 
odic boundary conditions includes not only single glue- 
ball states, but also states consisting of several glueballs 
and/or torelons (gluon excitations which wrap around 
the toroidal lattice), a means of identifying the single 
glueballs must be employed to prune away all of the other 
unwanted states. An additional small-volume simulation 
is done to assist in this identification and to study the 
systematic errors from finite volume. Finally, discretiza- 
tion errors are treated by extrapolating the energies to 
the continuum limit and determining the continuum spin 
quantum numbers. The end result is a nearly complete 
survey of the glueball spectrum in the pure gauge theory 
below 4 GeV. We find a total of thirteen glueballs; two 
other tentative candidates are also located. With the ex- 
ception of the light glueballs in the ++ sector, our results 
significantly improve upon those from previous studies 
of the complete low-lying glueball spectrum. 
This paper is organized as follows. The details of the 
simulations, including the construction of the glueball op- 
erators, the generation of the gauge-field configurations, 
the extraction of energies from Monte Carlo estimates of 
the correlation functions, and the lattice spacing determi- 
nations in terms of the hadronic scale vq , are described in 
Sec. ||. All of our energy estimates in terms of the inverse 
temporal lattice spacing are presented in this section. In 
Sec. |ni], the differentiation of single glueball states from 
two-glueball and torelon-pair states is discussed. System- 
atic errors from finite volume are studied in Sec. [V. The 



removal of lattice spacing errors, including the extrapo- 
lations to the continuum limit and the identification of 
the continuum spin quantum numbers, is described in 



i 



Sec. [v| We also discuss the problematical scalar states 
in this section and cite ongoing efforts to reduce their 
discretization errors. Sec. VI presents a discussion of the 



spectrum, and our findings are summarized with an out- 
line of future work in the concluding Sec. VII. 



II. SIMULATION DETAILS 

Our glueball mass determinations rely on numerical 
simulations of gluons on a hypercubic Euclidean space- 
time lattice with spatial and temporal spacings a s and 
a t , respectively. The gluons are described by the im- 
proved action Sn used in Ref. The couplings in the 
action depend on two parameters, f3 and £, and are de- 
termined using a combination of (tree-level) perturbation 
theory and mean-field theory, implemented by renormal- 
izing the spatial link variables Uj(x) — > Uj(x)/u s , where 
Us is given by the fourth root of the average plaquette 
|^| . The temporal link variables are not renormalized. 
The lattice anisotropy a s /at is given by £ at tree-level 
in perturbation theory. This action is intended for use 
with at <C a s , has 0(a A s ,a%,a s a 2 a ) discretization errors, 
where a s is the QCD coupling, and couples only nearest- 
neighbor time slices, ensuring the free-gluon propagator 
has no spurious modes. In all cases, glueball effective 
masses are seen to converge monotonically from above. 
This is a very desirable feature since it validates the use 
of variational techniques to minimize excited-state contri- 
butions to the glueball correlation functions. Such tech- 
niques are crucial for precise glueball mass extractions. 

On a simple cubic lattice, zero-momentum station- 
ary glue states are characterized by their transformation 
properties under the octahedral point group O, combined 
with parity and charge conjugation operations. O has 24 
elements (which are all proper rotations) that fall into 
five conjugacy classes; the single- valued irreducible repre- 
sentations are labeled A\, A2, E, Ti, and T2, (Schonfliess 
notation pi) and have dimensions 1, 1, 2, 3, and 3, re- 
spectively 71 . The inclusion of parity results in the sym- 
metry group known as Oh = O <g> C,, where Ci denotes 
the two-element group consisting of the identity opera- 
tion and spatial inversion. The conventional labels for 
the irreducible representations of Oh are obtained from 
those of O by appending a subscript g for representations 
corresponding to states which are even under parity and 
u for odd parity representations. However, we shall use 
a slightly different notation. Instead of the subscripts 
g and u, we use superscripts + and — , respectively, to 
indicate the eigenvalue P of parity. Glueball states are 
also eigenstates of charge conjugation. We denote the 
eigenvalue of charge-conjugation parity by C, as usual, 
and introduce an additional superscript in the represen- 
tation labels. We refer to the full symmetry group of 
zero-momentum glueball states on a simple cubic lattice 
as 0% or O pc ; the irreducible representations are labeled 



A PC , A PC , E PC , Tf c , and T 2 PC . For convenience, we 
use R when referring to these labels in general. Note that 
when we use one of these labels to identify a particular 
state, we refer to the lowest-lying zero-momentum state 
in the symmetry channel indicated by the representation 
label. The first-excited state in a particular symmetry 
channel will be denoted by the representation label with 
an asterisk. 

The mass of a glueball G having spin J, parity P, and 
charge-conjugation parity C can be extracted from the 
large-i behavior of a lattice-regulated correlation function 
C(t) = (0|$( fl )(*)<l (fl) (0)|0), where R is any irreducible 
representation of O^ occurring in the subduced represen- 
tation J PC lOg, and $<«)(*) = - (0|$( fl )(i)|0) 
is a gauge-invariant, translationally- invariant, vacuum- 
subtracted, real operator capable of creating a glueball 
from the QCD vacuum |0). As the temporal separa- 
tion t becomes large, this correlator tends to a single de- 
caying exponential limt—,00 C(t) — Z exp(— mat), where 
tuq is the energy of the lowest-lying state which can 
be created by the operator <§>( R \t). To determine mc, 
the correlator C(t) must be calculated for large enough 
t such that it is well approximated by its asymptotic 
form. Unfortunately, stochastic fluctuations in C(t) re- 
main roughly constant with t while the signal falls rapidly 
and hence, the use of a glueball operator for which C(t) 
attains its asymptotic form as quickly as possible is cru- 
cial for reliably extracting vtiq. The energies of ex- 
cited states in representation R can be obtained from 
the large-i behavior of a matrix of correlation functions 
Cij(t) = (0\^ R \t) ^ R) (0)|0), where each of the glueball 

operators &\ (t) transforms as R under all symmetry 
operations. Again, it is very important to use opera- 
tors for which the matrix elements CV, (t) attain their 
expected asymptotic forms for t as small as possible. 

Such operators can be constructed by exploiting link- 
smearing and variational techniques as previously de- 
scribed in Ref. [Q . For each irreducible representation R, 
glueball operators on a given time-slice are constructed 
in a sequence of three steps. First, a set of six smear- 
ing schemes are applied to the spatial link variables. 
Each scheme is a sequence of single-link and double-link 
mappings which depend on parameters X s and A/, re- 
spectively. We use the same six schemes described in 
Ref. . Secondly, a set of basic real operators <j>^ (t) is 
constructed using linear combinations of gauge-invariant, 
path-ordered products of the smeared link matrices about 
various closed spatial loops; each such linear combination 
is invariant under spatial translations and transforms ac- 
cording to the irreducible representation R. The loop- 
shapes employed in our calculation, shown in Fig. [j], are 
chosen for their ease of computation; all orientations of 
these operators can be computed very efficiently by first 
storing the untraced products of links around the twelve 
spatial plaquettes stemming from each site on a time- 
slice, then tracing the appropriate products of these ob- 
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TABLE I. The glueball simulation parameters. Values for 
the coupling /3, input aspect ratio parameter f , the mean-link 
parameter w S) the single- link smearing weight A s , the two- link 
smearing weight A/, and the lattice sizes are listed. Results 
for the hadronic scale ro in terms of the lattice spacing a s are 
also given. The approximate spatial lattice spacings a s are 
determined assuming r^ 1 — 410(20) MeV. 



£ u s A s A/ Lattice ro/a s 



1.7 5 0.295 0.1 0.5 (Fx 30 1.224(1) 

1.9 5 0.328 0.1 0.5 6 3 x 30 1.375(2) 

2.2 5 0.378 0.1 0.5 8 3 x 40 1.761(2) 

2.4 5 0.409 0.1 0.5 8 3 x 40 2.180(6) 

2.5 5 0.424 0.1 0.5 10 3 x 50 2.455(6) 



a s /ro a s (fm) 



0.8169(9) 0.39 

0.727(1) 0.35 

0.5680(5) 0.27 

0.459(1) 0.22 

0.407(1) 0.20 



3.0 3 0.500 0.4 0.5 15 3 x 45 4.130(24) 0.2421(14) 0.12 



FIG. 1. The Wilson loop shapes used in making the basic 
glueball operators. 



jccts. Both single and double windings around the paths 
are used; this allows us to double the number of raw oper- 
ators with only a small increase in computational effort. 
For each symmetry channel except , four irreducible 
combinations are then chosen and applied to the smeared 
links from each of the six schemes, yielding a total of 24 
basic operators in each channel. For the , only the 
last shape in Fig. [I] can be used and produces a total of 
12 basic operators. Lastly, the glueball operators & R \t) 
are formed from linear combinations of the basic oper- 
ators, $^(t) = ^„J)i (t), where the coefficients 

Va are determined using the variational method. This 
involves first obtaining Monte Carlo estimates of the large 
correlation matrix 



(1) 



where ^(t) = ^"'{t) - {O]^' (t)\0) . In practise, this 
vacuum subtraction is only performed for the A^ + chan- 
nel as the expectation value vanishes identically in all 



other sectors. The coefficients v. 
minimize the effective mass 



(R) 



are then chosen to 
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(2) 



where the time separation for optimization is fixed to 
t]j = 1. Other values of tjj are used as consistency 
checks. Let denote a column vector whose elements 



are the optimal values of the coefficients Va 
satisfies the eigenvalue equation 



(R) 



and the eigenvector Vq corresponding to the lowest ef- 

( R) 

fective mass mo(tn) yields the coefficients VqJ for the 
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This vector 



(3) 



operator $Q R \t) which, under ordinary circumstances, 
best overlaps the lowest-lying glueball Go in the channel 
of interest. A sequence of operators '(t), ^^(t), . . . 
which predominantly overlap excited glueball states can 
also be constructed using the higher-mass eigenvectors of 
Eq. (1). 

Monte Carlo estimates of the correlator matrix ele- 
ments given in Eq. ([!]) were obtained for all 20 irreducible 
representations in five simulations using an input aspect 
ratio parameter £ = 5. The values for the coupling /3, 
mean-link parameter u s , smearing weights A s and A / , and 
the lattice sizes used in these runs are listed in Table |. 
An additional run at a smaller lattice spacing (~ 0.12 fm) 
and using £ = 3 was done for the , E ++ , and T 2 t ~ + 
representations only. A smaller-a s measurement helped 
to obtain a reliable continuum-limit extrapolation for the 
troublesome state. The input parameters for this 

run are also given in Table |. All computations were 
carried out on DEC Alpha and Sun Ultrasparc work- 
stations. Configuration ensembles were generated us- 
ing Cabibbo-Marinari (CM) pseudo-heatbath and SU(2) 
subgroup over-relaxation (OR) methods. Link variables 
were updated in serial order on the lattice. Three com- 
pound sweeps were performed between measurements, 
where a compound sweep is one CM updating sweep fol- 
lowed by bor OR sweeps. The measurements were av- 
eraged into bins of ri m b, and n\>xas bins were obtained. 
For the (3 = 3.0, £ = 3 run, iior = 5, n m b = 50, and 
"bins = 80. For the (3 = 2.5, £ = 5 run, hor = 5, 
fi m b = 20, and n^ns = 318. For all of the other simu- 
lations, iior = 3, n m b = 100, and "bins = 100. Crude 
checks for residual autocorrelations were done by over- 
binning by factors of two and four; in all cases, statistical 
error estimates remained unchanged. 

In the final analysis phase, the glue energies vhq were 
extracted using a two-step procedure. First, the large 
correlation matrices in each channel were reduced to 
smaller 3x3 matrices Cab(S) for A, B = 0,1,2 using 



TABLE II. Glue energy estimates in terms ol a^ 1 for the five £ = 5 simula- 
tions. The levels are labeled by the irreducible representations of the cubic point 
group under which their corresponding stationary states transform. First-excited 
states in a representation are indicated by an asterisk. 





(3 = 1.7 


P = 1.9 


P = 2.2 


P = 2.4 


P = 2.5 




0.578(5) 


0.475(4) 


0.362(3) 


0.303(3) 


0.288(4) 


A*++ 


1.19(2) 


0.92(3) 


0.697(6) 


0.569(4) 


0.511(7) 


A ++ 

A 2 




1 971^ 








E++ 


0.924(8) 


0.844(6) 


0.667(4) 


0.538(3) 


0.472(4) 


E*++ 


1.29(2) 


1.093(12) 


0.878(7) 


0.723(9) 


0.652(8) 


T++ 


1.55(2) 


1.32(2) 


1.00(2) 


0.834(4) 


0.728(8) 


T*++ 


1.73(3) 


1.52(2) 


1.23(4) 


0.909(15) 


0.823(13) 


± 2 


1.103(8) 


0.918(7) 


0.686(4) 


0.542(2) 


0.477(3) 


rp* + + 
± 2 


1.41(2) 


1.228(12) 


0.938(5) 


0.730(8) 


0.660(6) 


A~ + 


1.31(2) 


1.06(5) 


756n A\ 




599('7 S ) 




1.86(9) 


1.47(4) 


1.08(2) 


0.836(9) 


0.72(2) 


A~ 2 + 





1.63(5) 


1.29(2) 


1.036(12) 


0.93(1) 


E-+ 


1.38(2) 


1.167(14) 


0.874(5) 


0.698(4) 


0.611(7) 


E*-+ 


1.85(5) 


1.49(3) 


1.085(9) 


0.890(6) 


0.782(13) 


T~ + 


1.56(2) 


1.42(2) 


1.155(9) 


0.873(13) 


0.82(1) 


T - + 

1 2 


1.297(14) 


1.148(10) 


0.882(5) 


0.695(4) 


0.619(3) 


rr + 


1.60(2) 


1.39(2) 


1.087(9) 


0.880(5) 


0.771(9) 






1.67(5) 


1.32(2) 


1.062(14) 


0.94(1) 


At~ 




1.29(2) 


0.999(11) 


0.796(7) 


0.700(14) 


E+- 




1.52(3) 


1.207(37) 


0.929(8) 


0.82(1) 


T+~ 


1.186(13) 


1.053(8) 


0.819(4) 


0.652(5) 


0.590(5) 


T*+~ 


1.55(3) 


1.297(14) 


1.025(8) 


0.794(9) 


0.73(2) 


A 2 




1.330(16) 


0.983(17) 


0.801(4) 


0.701(7) 


T~"H — 
1 2 




1.45(2) 


1.10(3) 


0.929(15) 


0.83(1) 


Ar 


1.83(7) 


1.61(5) 


1.354(20) 


1.04(3) 


0.98(1) 


A 2~ 




1.65(6) 


1.201(18) 


0.96(3) 


0.82(2) 


E— 


1.59(3) 


1.38(2) 


1.094(11) 


0.875(6) 


0.78(1) 


T~ 


1.72(4) 


1.46(2) 


1.07(3) 


0.877(5) 


0.760(9) 


T 2 - 


1.63(3) 


1.41(2) 


1.114(8) 


0.886(6) 


0.77(1) 



the variational coefficients of the three lowest mass eigen- 
states of Eq. (||) : 

C AB (t) = 5>|$<f>(r+*) ^(r)|0). (4) 

r 

Secondly, the expected large-t functional forms were fit 
to the elements of these optimized correlators. To ob- 
tain an estimate of the energy mc (in terms of a^ 1 ) 
of the lowest-lying state in each channel, a single expo- 
nential was fit to the ground-state correlator Coo(i) for 

i ^min; ■ • * ) ^max- 

Coo(t) = Z 00 {e-" iG o* + e - mG o(^-*)} , ( 5 ) 

where T was the temporal extent of the periodic lattice. 
To determine the energies tog p of the excited states and 



another estimate of mc , the M xM optimized correlator 
matrix was also fit for t = t min , . . . , i max using the form 

A/-1 

Cab® = Za p Zb p { e ~ r " Gpt + e"^^} , (6) 

p=0 

for M — 2,3. Various fit regions t m - ln to i max were used 
to check for consistency in the extracted values for the 
masses. Best-fit values were obtained using the corre- 
lated x 2 method. Error estimates were calculated using 
a 1024— point bootstrap procedure; in all cases, error es- 
timates were nearly symmetric about the central best-fit 
values and were thus averaged to simplify presentation. 
Our fits are far too numerous to list here; additional de- 
tails are available from the authors upon request. 

Our final estimates of the glue energies in terms of a^ 1 
for the £ = 5 simulations are presented in Table [□]; energy 
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TABLE III. Glue energy estimates in units of a t 1 from 
the f3 = 3.0, £ = 3 simulation. 





0.318(4) 


E++ 


0.476(3) 


T++ 


0.476(3) 



estimates from the £ = 3 run are listed in Table II] . In or- 
der to convert these results into physical units, the lattice 
spacings at must be determined for each simulation. The 
hadronic scale parameter ro § defined in terms of the 
force between static quarks by [r 2 dV(r)/dr] r=ra — 1.65, 
where V(r) is the static-quark potential, is a useful quan- 
tity for this purpose. The values for ro in terms of o s cor- 
responding to each glueball simulation were determined 
by measuring the static-quark potential in separate sim- 
ulations. The results are listed in Table |. Further details 
concerning the calculation of ro ja s are given in Ref. |Jp . 
Note that in computing ro/a s , the input value £ was used 
for the aspect ratio a s j at - The consequences of doing so 



are discussed in Sec. VA 



III. GLUEBALL IDENTIFICATION 

The spectrum in a box with periodic boundary condi- 
tions includes not only single glueball states, but also 
states consisting of several glueballs and/or torelons 
(gluon excitations which wrap around the toroidal lat- 
tice). We expect that the operators used in our cor- 
relators will couple most strongly to the single glueball 
states, but we cannot be certain that mixings with the 
multi-glueball and torelon states will be negligible. Re- 
call that the asymptotic behavior of the correlation func- 
tion associated with an operator is dominated by 
the lowest-lying eigenstate which mixes with $. If a 
multi-glueball or torelon state has a lower energy than 
the lightest glueball in a given symmetry channel, then 
the possibility exists that the energy we extract from the 
asymptotic decay of the correlator will be that of the 
multi-glueball or torelon state. Thus, a means of differ- 
entiating the single glueball states from all other states 
is required. 

First, given mass estimates of the lowest-lying few glue- 
balls, one can easily determine the approximate locations 
in a given symmetry channel of the two-glueball states 
having zero total momentum. If the simulation results in 
that channel lie significantly below the two-glueball en- 
ergy estimates, one can almost certainly rule out a multi- 
glueball interpretation. 

Secondly, one can study the manner in which each en- 
ergy level changes as the lattice volume is varied. The 
energy of a single glueball state depends on the lattice 



volume in a markedly different way from that of a multi- 
glueball or torelon state. 

A third possibility is to include additional operators 
in the correlation matrices which are expected to couple 
much more strongly with the two-glueball and torelon 
states. The construction of operators which best overlap 
the lowest-lying eigenstates of interest using the varia- 
tional method then involves not only single glueball op- 
erators, but also the new two-glueball and torelon op- 
erators. The coefficients obtained from the variational 
optimization can be used to estimate the mixings of the 
additional operators with the low-lying eigenstates of in- 
terest. If the mixings of the additional two-glueball and 
torelon operators with an eigenstate are very small, a 
single glueball interpretation is assured; in such a case, 
the addition of the new operators does not affect the 
extracted energy. If the mixings are significant, the addi- 
tion of the two-glueball/torelon operators will lower the 
extracted energy, ruling out a single glueball interpreta- 
tion. 

Ideally, it would be best to apply all three of these 
methods. However, for this initial scan of the glueball 
spectrum, we decided for reasons of simplicity to rely 
mainly on the first method. Having obtained the lowest- 
lying one or two states in each symmetry channel, we 
determined the approximate locations of the two-glueball 
and torelon states. Simulation results lying near or above 
these thresholds were then excluded from further consid- 
eration. In other words, we used the two-glueball and 
torelon thresholds as filters to remove possibly extrane- 
ous states. One additional simulation was done to study 
systematic errors from finite volume. This run also served 
to confirm the single glueball nature of the states lying 
well below the two-glueball/torelon thresholds. Note that 
the levels lying near or above these thresholds cannot be 
ruled out as single glueball states; rather, one can say 
only that the interpretation of such states requires addi- 
tional information. 



A. Two-glueball states 

In order to identify the genuine single glueball states, 
we first determined the approximate locations of the two- 
particle states using the mass estimates of the lightest few 
glueballs. In estimating these locations, we assumed that 
the energy of a two-glueball state was given by 



p 1 



P 2 



(7) 



where mi is the rest mass of the glueball having mo- 
mentum p and mi is the rest mass of the other glueball 
which has momentum —p. Note that on a periodic lat- 
tice having N s sites in each of the three spatial directions, 
the allowed momenta are discrete p = 2Tr(n Xy n y ,n z )/L, 



where L = a*N* and 



and n z are integers satis- 



fying —N s /2 < rij < N s /2. The above energy estimates 



TABLE IV. The little groups of p= 2ir(n x , n y , n z )/L on a 
simple cubic lattice with periodic boundary conditions. Note 
that l,m,n/0 and I 7^ m, m 7^ n, and I 7^ n. 



(n x ,n y ,n z ) Little group 

(0,0,0) o~ h 

(0,0, n) C iv 

(0,n, n) Civ 

(n,n,n) Csu 

(0,m, n) C s 

(m,m,n) C s 

(l,m,n) Ci 



neglect all interactions between the two glueballs; this 
should not introduce serious error since we expect these 
interactions to be short-ranged, being mediated by scalar 
glueball exchange at large distances. Eq. (^) also assumes 
that the rest masses and dispersion relations of the propa- 
gating glueballs are unaffected by finite volume and finite 
lattice spacing errors. We have verified this assumption 
in the case of the scalar glueball for n^+riy +n 2 z < 9 on 
an 8 3 x 40 lattice at (3 = 2.4 and £ = 5. 

To facilitate our discussion of the two-glueball states, 
we first point out some features of single, propagating 
glueballs in a finite box with periodic boundary condi- 
tions. Here, glueball states are characterized by their 
transformation properties under 0\ c , the simple cu- 
bic crystallographic space group 0\ extended to include 
charge-conjugation. The group 0\ is isomorphic to the 
semi-direct product of the group of pure (discrete) trans- 
lations and the group of pure (discrete) rotations and 
reflections about a given center. Thus, a propagating 
glueball state may be classified by its momentum p and 
by its transformation properties under the subgroup of 
0\ which leaves p invariant (the little group of p) . The 
little groups corresponding to various momentum orien- 
tations are listed in Table |v|. Hence, the irreducible 
representations of the little group may be used to iden- 
tify propagating glueball states. The little group varies 
with the momentum orientation. This means that the 
partitioning of the physical glueball states into the irre- 
ducible representations of the lattice symmetry groups 
differs depending on the momentum orientation. For ex- 
ample, consider the 2 ++ glueball. When at rest, three of 
the five polarizations of this glueball appear in the 
representation of 09, the little group of p = (0, 0, 0), and 
two of its polarizations occur in the E ++ representation. 
When p = (0, 0,p) for p^O, the five polarizations of the 
2 ++ glueball (which are no longer eigenstates of parity) 
split across the , B± , £?^~ , and E + representations of 
the little group C% v . 

We determined the lowest-lying two-glueball states in 
each symmetry channel by repeating the following se- 
quence of steps for each allowed momentum vector p and 



each choice of two glueballs G\ and Gi- For the mo- 
ment, assume that G\ and G2 are distinguishable. Note 
that G\ and Gi refer to the irreducible representations 
of O c , the little group for zero momentum. First, we 
identified the little group L(p) of p. Secondly, the char- 
acters ^( Gl ) and ^( G2 ) of the representations G\ and 
G2 were subduced into the little group, yielding the 
characters x^ G ^ [ L(p), which were then decomposed 
into the irreducible representations fj, of the little group: 
yAGj) 1 L(p) = cjti Next, for each [i\ such 

(G ) (G ) 

that a;,/ 7^ and each ^ 2 such that c/t 2 7^ 0, we 
formed the direct product ^(aw®^) ^ obtain the char- 
acter corresponding to the two-glueball state. Since the 
total momentum of this two-glueball state is zero, it can 
be characterized by the irreducible representations of O c . 
A representation of O c was formed by constructing a set 
of coset representatives and applying the method of in- 
duction, and the induced character ^(mi®a«2)-|-qC wag g_ 
nally decomposed into the irreducible representations of 
O c . When G\ and G2 were indistinguishable, the above 
procedure was modified to include Bose symmetrization. 

This procedure was carried out using the rest energies 
obtained in the (3 = 2.5, £ — 5 simulation. The lowest- 
lying two-glueball states in each symmetry channel are 
listed in Table [v|. These levels, along with all higher 
lying levels, are shown as dashed line segments in the 
shaded region in Fig. ^. Any energy lying well below the 
shaded region in this figure can be safely interpreted as 
a single glueball state (these are indicated by the black- 
outlined green-filled boxes). States lying slightly below 
(the orange boxes with no outlines) or above (the red 
boxes) the two-glueball thresholds must be regarded with 
caution. Again, we remind the reader that we cannot 
rule out a single glueball interpretation for these levels; 
additional information is needed to determine the nature 
of these states. Since it is not our intent in this paper 
to obtain such information, we exclude these levels from 
the spectrum of single glueball states for the time being. 

B. Torelon pairs 

Torclons are gluonic excitations which wind around the 
periodic boundaries of the lattice. They may be classi- 
fied according to their behavior under a set of discrete Z% 
symmetries of the SU(3) pure gauge theory. The gauge 
action is invariant under multiplication of every link in 
the [i direction (/i = 0, 1, 2, 3) at a given /x coordinate by 
the same member of Z3, the center of SU(3). The torelon 
is an eigenstate of the transfer matrix which transforms 
nontrivially under such symmetry operations. For exam- 
ple, the spectrum of glue on a periodic lattice contains 
a torelon eigenstate which picks up a phase exp(27n/3) 
and a state which picks up a phase exp(— 27ri/3) under 
the multiplication of every link in the x-direction at a 
given x coordinate by the center member exp(27rz/3). In 
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TABLE V. The lowest lying state consisting of two free 
glueballs in each symmetry channel for the /3 = 2.5, £ = 5 
simulation. Each state is comprised of a glueball having a 
momentum p = 2ir(n x ,n y ,n z )/L, where L = 10a 3 , and an- 
other glueball having momentum —p. 



Channel 


Glueballs 


(n^, rty, 712) 




(At + ,At + ) 


(0, 0, 0) 


At + 


(A++,E + +) 


(0, 0, 1) 


E++ 


(At\At + ) 


(0, 0, 1) 


T++ 


(At + ,A++) 


(0, 1, 2) 


T++ 


( A++ A++\ 

(A . A 1 ) 


(0, 1, 1) 




(Al + ,Af+) 


(0, 0, 0) 


A- 2 + 


(A++.T++) 


(0, 0, 1) 


E-+ 


(A++.T++) 


(0, 0, 1) 


T - + 


(A++,i5++) 


(0, 0, 1) 




(A++,i5++) 


(0, 0, 1) 


At~ 


(A++.T+-) 


(0, 1, 2) 


At~ 


(A++.T+-) 


(0, 1, 1) 


E+- 


(A++.T+-) 


(0, 1, 1) 


T+- 


(A++.T+-) 


(0, 0, 0) 


T+- 


(A++.T+-) 


(0, 0, 1) 


Ar 


(A++.T+-) 


(0, 0, 1) 


Ar 




(0, 1, 1) 


E— 


(A++.T+-) 


(0, 0, 1) 






(0, 0, 1) 


T 2 — 




(0, 0, 1) 



fact, there are three such pairs of modes corresponding 
to the x, y, and z directions. These torelons can also 
have a center-of-mass momentum in the two spatial di- 
rections transverse to their flux direction. For large L, 
a torelon at rest has an energy given approximately by 
aL, where a is the string tension from the static-quark 
potential and L is the spatial extent of the lattice. We 
have confirmed this in a simulation on an 8 3 x 24 lattice 
at (3 = 2.4, £ = 3. Hence, the torelon mass is strongly 
dependent on the volume of the lattice. 

Our glueball operators, being closed Wilson loops 
which do not wrap around the lattice, are invariant un- 
der these center symmetry transformations. This means 
that the glueball operators cannot create a single torelon 
state, but the creation of torelon pairs of opposite center 
charge is possible. If the two torelons of opposite charge 
do not substantially interact, then the lowest energy of 
such a state is 2aL, which has the value 0.9 when placed 
on Fig. ||. This lies sufficiently high to discard from con- 
sideration, even for the A~{ state since a state of two 
torelons of opposite center charge and total zero mo- 
mentum must be symmetric under charge conjugation. 
Fortunately, a torelon pair can be easily detected in a 
finite-volume study since their energy depends strongly 
on L. An additional simulation was done to measure the 
changes in all energy levels as the lattice volume was re- 
duced. The results of this simulation are presented in the 
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FIG. 2. Comparison of the pure-glue spectrum obtained 
from the (3 = 2.5, £ = 5 simulation to the approximate loca- 
tions of the two glueball states. The boxes are the simulation 
results; the standard deviations in these mass estimates are 
indicated by the vertical heights of the boxes. The dotted 
line segments shown in the upper shaded region indicate the 
approximate locations of states consisting of two free glue- 
balls having zero total momentum. All energies are in terms 
of Of . The representations of the cubic point group which 
label the states are indicated along the horizontal axis. The 
most likely J PC interpretations of the states are also shown. 

next section. No energy reductions of sufficient magni- 
tude were found to suggest that any of our states could 
be interpreted as a torelon pair. 

IV. FINITE VOLUME ERRORS 

An additional simulation was done to measure the sys- 
tematic errors in our results from finite volume: a run at 
(3 = 2.4, £ = 5 on a 6 3 x 40 lattice. The original 8 3 x 40 
lattice has a spatial volume of (1.76 fm) 3 , assuming that 
a s ~ 0.22 fm from r^ 1 = 410(20) MeV. The additional 
simulation measures the changes in the glueball masses 
as the volume is reduced from (1.76 fm) 3 to (1.32 fm) 3 . 
The input parameters used in this small- volume run were 
the same as those used in the larger volume simulation. 

Let niQ denote the energy of a state G on the origi- 
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FIG. 3. Finite-volume effects on the results of the (3 = 2.4, 
£ = 5 simulation. Each point shows the fractional change 
5g = 1 — rn G /mG in the energy of a stationary state G, where 
tug is the energy of G as measured on an 8 3 x 40 lattice, and 
m,Q is the energy of G as measured on a smaller 6 3 x 40 lattice. 
The state G corresponding to a given point is specified by 
combining the representation label below the point along the 
horizontal axis with the PC value shown to its left along the 
vertical axis. The solid lines indicate Sg = 0, the dotted lines 
above the solid lines indicate Sg = 0.02, and the dotted lines 
below the solid lines indicate Sg = —0.02. 



nal 8 3 x 40 lattice, and m G denote the energy of G as 
measured on the smaller 6 3 x 40 lattice. The fractional 
change in the energy is defined by Sg — l — m G /rriG- The 
results for these fractional changes are shown in Fig. ||. 
Each point shows the fractional change in the energy of 
a state G specified by combining the representation la- 
bel below the point on the horizontal axis with the PC 
value shown to its left on the vertical axis. The solid lines 
indicate Sg = 0, the dotted lines above the solid lines in- 
dicate Sq — 0.02, and the dotted lines below the solid 
lines indicate Sg = —0.02. The largest effects from finite 
volume occur in the A* ++ and T* + ~ states. All other 
changes are statistically consistent with zero, indicating 
that systematic errors in these results from finite vol- 
ume are negligible. These results confirm the single glue- 
ball nature of all of the states, with the exception of the 
AJ ++ , lying well below the two-glueball thresholds (the 
black-outlined green-filled boxes in Fig. ^). Although the 
change in the energy of the A 1 ++ is not very large, it is 
sufficient to warrant further study of this level. For this 



reason, we withhold judgment on whether or not this 
level is a single glueball. Note that most of the states ly- 
ing near or above the two-glueball thresholds (the orange 
and red boxes in Fig. ||) show very little finite volume de- 
pendence, suggesting that these states might actually be 
long-lived glueball resonances. Further study would be 
required to resolve this issue. 



V. LATTICE SPACING ERRORS 

There are two aspects to the removal of systematic 
errors from finite lattice spacing: extrapolation of the re- 
sults to the continuum limit a s — ► O, and the identification 
of the continuum spin quantum numbers. In this section, 
we first carry out the a s — > extrapolations of the candi- 
date single-glueball levels remaining after the analysis of 
Sec. Ill, then deduce their continuum spin content. 



A. Continuum limit extrapolations 

The glueball mass estimates in terms of a^" 1 , listed in 
Table ||, were combined with the determinations of the 
hadronic scale r$/a s presented in Table [|. The results 
are shown in Figs. (4||7|. In these figures, the dimension- 
less products of ro and the glueball mass estimates are 
shown as functions of (a s /ro) 2 . The solid symbols indi- 
cate the results from the £ = 5 simulations. The open 
symbols appearing to the right of the vertical dashed lines 
indicate results from the (3 = 3.0, £ = 3 run, as well as 
£ = 3 results for the A~l + and A* ++ channels previously 
obtained in Ref . |ffl] . To remove discretization errors from 
our glueball mass estimates, the results for each level in 
these figures must be extrapolated to the continuum limit 
a s /r — >0. The discretization errors can then be seen as 
the deviations of the finite-a s results from these limiting 
values. 

From perturbation theory, the leading discretization 
errors in our action are expected to be 0(af,a^,a s a,s). 
The agreement of the A++ , A*++, E++, T++\ and T+~ 
glueball masses extracted using £ = 3 and £ = 5 (see 
Figs. 12 and 14 in Ref. Q) suggests that the 0(af) er- 
rors are negligible. Some evidence for the smallness of 
the 0(a s a 2 s ) errors comes from our earlier glueball sim- 
ulations H which used the one-loop improved Liischer- 
Weisz action |l(J. After tadpole improvement, the ra- 
diative corrections to the couplings in this action were 
found to be very small. For these reasons, we expect 
that the 0(a^,a s a 2 ) errors will be negligible compared 
to the O(cig) errors. 

In assuming a s /a t = £ (where £ is the input anisotropy 
parameter in the action), we introduce 0(a s ) errors in 
our estimates of the glueball masses multiplied by ro- 
(Note that these errors do not enter into ratios of the 
glueball masses.) Such errors can be avoided by instead 
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FIG. 4. Mass estimates of the PC = ++ glueballs in 
terms of ro against the lattice spacing (a a /ro) 2 . The solid 
symbols indicate results from the £ = 5 simulations, and the 
open symbols on the right-hand side of the vertical dashed 
line indicate results from the £ = 3 simulations. The solid 
curves are best fits to the simulation results for each state 
using ipi (a s ) from Eq. ( Jic| ) for the Af + and A^ ++ levels and 
tpo (a s ) from Eq. (^j) for all other levels. The open symbols on 
the left-hand side of the vertical dashed line show the extrap- 
olations to the continuum limit using the best-fit forms. 

setting a s /at in a more physically motivated fashion, such 
as by comparing the spatial and temporal length scales 
extracted from appropriate correlation functions. We can 
also use perturbation theory to adjust the couplings in 
our action to remove these errors order by order in a s . 
However, we estimated the errors caused by imposing 
a s /at — £ (see below) and found them to be too small 
to warrant the additional complexity of another a s /at- 
setting scheme or the effort required to calculate the one- 
loop corrections to the action. A simpler approach is to 
incorporate the 0(a s ) anisotropy errors into our contin- 
uum limit extrapolations. Unfortunately, their depen- 
dence on a s is not well known and, as we shall see, they 
are generally smaller than the much more rapidly vary- 
ing 0(a A s ) errors. Detecting their effects in a fit to about 
five data points for a s from 0.2 to 0.4 fm is not feasi- 
ble. Thus, we decided to adopt the following approach: 
to extrapolate assuming 0{a i s ) errors only, then include 
a systematic uncertainty in our continuum-limit results 
from the 0(a s ) anisotropy errors. 
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FIG. 5. Mass estimates (solid symbols) of the PC = — + 
glueballs in terms of ro against the lattice spacing (a s /ro) 2 - 
The solid curves are best fits of (po(a s ) from Eq. M to the 
results for each state; the open symbols are the continuum 
limit extrapolations. 

One way to estimate this uncertainty is to compare 
measurements of the static-quark potential extracted 
from Wilson loops taken along the different spatial and 
temporal axes of the lattice [1TJ . The anisotropy errors 
can be quantified by defining a s /at in terms of these 
different potentials and comparing the result to £. If 
we denote the determination of the aspect ratio from 
the potentials by [a s /at]v and define by the rela- 
tion [a s /ot]y = Z^, then the deviation of Z^ from unity 
gives us a measure of the fractional error from assum- 
ing a s /at = £• The effect of these errors is to modify 
the multiplicative ro/dt factors used to convert the sim- 
ulation results given in terms of a^ 1 into units of r^ 1 
suitable for extrapolation. Using the functional depen- 
dence of ro I at on £ from a fit to the static-quark poten- 
tial, we determine that (Z^ — l)/2 gives us an estimate 
of the fractional uncertainty in our continuum limit re- 
sults from the aspect ratio errors. Without mean-link 
improvement, Z^ can deviate from unity by as much as 
30%. When the action includes mean-link factors, the 
corrections are found to be small, typically a few per 
cent. For example, for the j3 — 2.4, £ = 5 run, we ob- 
tained Z^ = 0.987(8); for (3 — 3.0, £ = 3, an estimate 
of Zc = 0.99(1) was found. If we assign a conservative 
2% error from Z%, then this amounts to a 1% system- 
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FIG. 6. Mass estimates (solid symbols) of the PC = + — 
glueballs in terms of ro against the lattice spacing (a s /ro) 2 - 
The solid curves are best fits of (po(a s ) from Eq. (bj) to the 
results for each state; the open symbols are the continuum 
limit extrapolations. 
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FIG. 7. Mass estimates (solid symbols) of the PC = 

glueballs in terms of ro against the lattice spacing (a s /ro) 2 - 
The solid curves are best fits of (po(a s ) from Eq. (fej) to the 
results for each state; the open symbols are the continuum 
limit extrapolations. 



atic uncertainty in our continuum-limit results from the 
anisotropy errors. 

Another way to estimate the errors due to the as- 
pect ratio is from the perturbative calculation of us- 
ing the anisotropic Wilson action (the analogous cal- 
culation using Sn has not yet been done). Modify- 
ing the results from Ref. [[L2| to include tadpole im- 
provement factors and writing the mean-link parameter 
u s = 1 — ev s ui 2 ^(£) + 0(a 2 ), one finds 



1 



{2tt MO-MO] + «< 2) (0} +0(a 2 s ), (8) 



where the values for c T (£) and c CT (£) can be obtained from 
Fig. 1 in Ref. From this equation, one sees that 

when £ = 3 — 5 and a s ~ 0.2, the aspect ratio receives 
less than a 2% correction; if the tadpole improvement 
factor u s (0 is omitted, a large 30% correction is found. 
For the improved action Sn, we expect that these values 
should be somewhat smaller. Again, we can assign a 
conservative 2% error in Z^ to obtain a 1% systematic 
uncertainty in our continuum estimates. 

To summarize, our approach is to extrapolate to the 
continuum limit using 



ip (a s ) = r m G + c 4 -4, 



(9) 



where r^mc and C4 are the best-fit parameters, then add 
a 1% systematic uncertainty from the 0{a s ) anisotropy 
errors. Eq. (p|) worked well in all cases except for the 
,4^+ and Ay*~ levels. The best-fit curves using Eq. @ 
are shown in Figs. the extrapolations of these curves 
to the continuum limit arc indicated in these figures by 
the open symbols on the left-hand sides of the vertical 
dashed lines. Note that the extrapolation uncertainties 
shown in these plots do not yet include the systematic 
anisotropy error. 

As discussed in Ref. [0, our results for the Af + 
and levels remain problematical. These levels 

have large finite-lattice-spacing errors which do not obey 
Eq. (||). There is growing evidence that these large dis- 
cretization errors are due to the presence of a critical 
end point of a line of phase transitions (not correspond- 
ing to any physical transition found in QCD) in the 
fundamental-adjoint coupling plane Jl3|-p^[. It has been 
conjectured that this critical point defines the continuum 
limit of a <p 4 scalar field theory |0. As one nears this 
critical point, the coherence length in the scalar channel 
becomes large, which means that the mass gap in this 
channel becomes small; all other observables, including 
glueball masses, appear to be affected to a much lesser 
extent. The effect of this critical point on the functional 



1 n 



form of the discretization errors in the scalar glueball 
mass is not well known, so we must proceed somewhat 
empirically. One possibility is that the critical point 
might enhance the perturbative 0(a s ) errors; another 
is that it may induce an 0{a 1 s ) error which would not 
otherwise be present. We fit several simple functions to 
the mass results and found that the following two 
functions worked very well: 



ipi(a s ) = r m G + c 4 -f 



-[d2^ + d i ^)[Ha s A)]- 



i \ a s 

f2{a s ) = r m G + c 2 -j 



a s a s 
C4 —r + C 6 -g , 



(10) 
(11) 



where romp, C4, Ce, cfe, and g?4 are the best-fit parame- 
ters. Eq. ( jll| ) is simply a cubic polynomial in (a s /ro) 2 . 
Eq. ( |l0| ) incorporates the expected leading dependence 
on the QCD coupling a s {a s ) ~ — 1/ ln(a s A) up to 0(a A s ). 
Various estimates of the QCD scale parameter A^g sug- 
gest that r A 
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0.6. Hence, we used r A = 0.5 and 
verified that our continuum limit estimates were insensi- 
tive to the choice of r$k in the range from about 0.3 to 
0.8. Both the di and d± terms were needed to achieve 
this insensitivity. 

The best-fit curves using Eq. (]l(]) are shown in Fig. ||. 
The best fit to the A+ + results has xV^df = 0.57, 
and for the we find x 2 /^ DF = 0.35. The ex- 

trapolation of these curves to the continuum limit yields 
r m{A+ + ) = 4.21(7) and r Q m{A* 1 ++ ) = 6.50(44). Using 
Eq. @), we obtain r m(A+ + ) = 4.30(8) with xV^df = 
0.53 and r m{A*+ + ) = 6.52(54) with xV^df = 0.35 
Since f\{a s ) was more closely connected to a pertur- 
bative analysis, we chose the estimates obtained using 
Eq. ( |io| ) for our final results, but added the differences be- 
tween the two extrapolations as a systematic error. After 
including an additional 1% anisotropy error, we end up 
with r Q m(A+ + ) = 4.21(12) and r m(A* ++ ) = 6.50(45). 

Note that the A^ + estimate differs slightly from our 
earlier estimate of 3.98(15) given in Rcf. [Q . Our previous 
extrapolation suffered from the absence of a mass mea- 
surement at a lattice spacing smaller than 0.2 fm; the 
need for such a measurement to obtain a reliable con- 
tinuum limit estimate for this level was acknowledged 
in Ref. [lj. The current study includes a new measure- 
ment at a lattice spacing near 0.1 fm; the inclusion of 
this new measurement is responsible for the slight differ- 
ence in the two extrapolations. Note that our improved 
estimate 4.21(12) agrees very well with the value 4.33(5) 
obtained by extrapolating existing Wilson action data for 
the scalar glueball mass. Unfortunately, the mass of the 
is very poorly determined because a measurement 
of the mass was not obtained in the ft = 3.0, £ = 3 

simulation. 

Recently, we have demonstrated that the discretization 



TABLE VI. Number of times each irreducible representa- 
tion of the octahedral group O occurs in the subduced rep- 
resentations J I O of the rotation group SO(3) restricted to 
subgroup O. 
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errors in the scalar glueball mass can be dramatically re- 
duced by simulating with an action which includes an 
additional two-plaquette interaction Jl5| |. With such an 
action, we should be able to substantially improve upon 
our determinations of the scalar-channel glueball masses 
in the near future. A study in SU(2) lattice gauge the- 
ory |l(| has also shown that lattice-spacing errors in the 
scalar glueball can be reduced by using the mean spa- 
tial and temporal links in Landau gauge for the values of 
the link variable renormalization parameters u s and ut, 
respectively. 



B. Spin identification 

The last step in our calculation of the glueball spec- 
trum is to identify the continuum spin content of each 
level. This is done by matching the observed patterns of 
degeneracies in the levels from different representa- 
tions to those expected for the various continuum J PC 
states. For example, a J = state occurs only in the Ai 
representation of O, a J = 2 state occurs in both the E 
and T 2 representations, and a state of J = 3 gets split 
across the A2, Ti, and T 2 representations. The numbers 
of times that the irreducible representations of the oc- 
tahedral group O occur in the subduced representations 
J I O of the rotation group SO (3) restricted to the sub- 
Given the values in this 



group O are listed in Table |VI , 
table and either the continuum limit estimates in Figs. 
or the results shown in Fig. [|, we can then deduce the 
continuum spin quantum numbers. 

Consider first the PC = ++ sector. The A± + state 
is not degenerate with any other level; hence, it can be 
identified as a J — state. The E ++ and T 2 ++ states 
are degenerate, implying that they correspond to the five 
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polarizations of a J = 2 glueball. The state is 

seen to have no degenerate partners, suggesting a J = 
excited state. The + state can correspond to J = 
3, 6, 7, 9, . . .. For all of these J values, there should be an 
accompanying level in the + channel, and such a level 
is observed. We conclude that the and T 1 ++ states 
correspond most likely to a J — 3 state, but less likely 
J = 6, 7, 9, . . . interpretations cannot be ruled out. The 
J = 3 assignment is also supported by model predictions, 



discussed in Sec. |VI 
In the PC 



+ sector, the A l + and A\ + states are 



states, and the degenerate 
2 glue- 

+ n „A T* — I- 



T 2 + states must correspond to a J 



easily identified with J 

E~+ and rr ~ + 

ball. The degenerate E* h and T 2 * h states most likely 
correspond to a J = 2 state as well, although as noted in 
Sec. Ill, the proximity of these levels to the two-glueball 
threshold leaves their status somewhat uncertain. Also, 
we cannot statistically rule out the possibility that, in 
combination with two states in the T 1 _+ irreducible rep- 
resentation, they are associated with a J — 5 glueball. 
We cannot rule out an accompanying degenerate level 
in the channel. If such a state exists, then these 

could be the levels corresponding toaJ = 4oraJ = 8 
glueball. 

The state must correspond to a J = 1 glueball, 
and the degenerate A\ , , and T£ levels corre- 

spond to a J = 3 glueball. A glueball having J = 2 is 
the most likely interpretation for the degenerate E~* and 
Tjf states. However, we cannot rule out the possibility 
of accompanying levels in the T± and A\ channels. 
Taking all possibilities into account, the alternate inter- 
pretations are J = 5, 7, 11. The very high lying A\ can 
be interpreted as a J = glueball, but J = 4, 6,8, . . . 
cannot be excluded. 

Finally, consider the PC = sector. The most prob- 
able scenario is as follows: the T-j" corresponds to a 
J = 1 glueball, the degenerate E~~ and T 2 ~~ are the 
five polarizations of a J = 2 glueball, and the A^ is 
a J = 3 state. Of course, we cannot rule out the pre- 
sumably higher-lying J = 6,7,9,11,... interpretations 
for the A 2 . Another possibility is that the E~~, T 1 , 
and T 2 are degenerate, in which case they could cor- 
respond to a J = 5 glueball. Even less likely is that 
all four levels are degenerate. In this case, one could 
interpret them as a single J = 7orJ=ll state, or 
as accidentally degenerate J = 3 and J = 2 glueballs. 
However, all higher-spin interpretations would require a 
degenerate ground state in cither the T" or T 2 ~~ chan- 
nel. Our correlator fits suggest these degeneracies are 
missing, making J = 5, 6, 7 . . . interpretations unlikely. 



VI. RESULTS AND DISCUSSION 

Our final results for the glueball spectrum in terms 
of r are given in Table VII. In this table, we assumed 



TABLE VII. Final continuum-limit glueball mass esti- 
mates ma- When a unique J interpretation for a state cannot 
be made, the other possibilities are indicated in the second 
column. States whose interpretation requires further study 
are indicated by a dagger. In column 3, the first error is 
the statistical uncertainty from the continuum-limit extrapo- 
lation and the second is the estimated uncertainty from the 
anisotropy. In the final column, the first error comes from the 
combined uncertainties in roma, the second from the uncer- 
tainty in 7-q 1 = 410(20) MeV. 





Other J 


rorriG 




m G (MeV) 


0++ 




4.21 (11) 


(4) 


1730 (50) 


(80) 


2++ 




5.85 (2) 


(6) 


2400 (25) 


(120) 


o-+ 




6.33 (7) 


(6) 


2590 (40) 


(130) 






6.50 (44) 


(7)t 


2670 (180) 


(130) 


1 + - 




7.18 (4) 


(7) 


2940 (30) 


(140) 


2"+ 




7.55 (3) 


(8) 


3100 (30) 


(150) 


3+" 




8.66 (4) 


(9) 


3550 (40) 


(170) 


*-+ 




8.88 (11) 


(9) 


3640 (60) 


(180) 


3++ 


6,7,9,... 


8.99 (4) 


(9) 


3690 (40) 


(180) 


\— 


3,5,7,... 


9.40 (6) 


(9) 


3850 (50) 


(190) 


2 *-+ 


4,5,8,... 


9.50 (4) 


(9)t 


3890 (40) 


(190) 


2"" 


3,5,7,... 


9.59 (4) 


(10) 


3930 (40) 


(190) 


3~ 


6,7,9,... 


10.06 (21) 


(10) 


4130 (90) 


(200) 


2+" 


5,7,11,... 


10.10 (7) 


(10) 


4140 (50) 


(200) 


o+- 


4,6,8,... 


11.57 (12) 


(12) 


4740 (70) 


(230) 



the most likely spin interpretations as described in the 
previous section and accordingly combined the contin- 
uum limit extrapolations (shown as open symbols on the 
left-hand sides of the vertical dashed lines in Figs. f|- 
R), then added the 1% anisotropy error, to obtain final 
estimates for the glueball masses in terms of ro. The 
combinations used to obtain these estimates are also 
indicated in Figs. by the J PC labels near the left 
vertical axes. Whe rever applicable, we have also indi- 
cated in Table VII any alternative spin interpretations 
which cannot be ruled out. These final estimates are also 
shown in Fig. ||. The 0* ++ and 2* ^ states are shown 



as dashed hollow boxes to indicate that their interpreta- 
tions as glueballs are tentative. Our concern about the 
0* ++ state stems from its non-negligible finite volume ef- 
fects; for the 2* h state, its nearness to the two-glueball 
threshold in our simulations is worrisome. Note that our 
estimates r o m(0++) = 4.21(12) and r m(2++) = 5.85(6) 
agree very well with 4.33(5) and 6.0(1), respectively, ob- 
tained by extrapolating the Wilson action simulation re- 
sults from Refs. |^,||,^7|jl8f| to the continuum limit. Sev- 



eral glueball mass ratios are presented in Table VIII. We 
can determine these ratios very accurately since, as noted 
earlier, they are not contaminate d by anisotropy errors. 
The uncertainties given in Table VIIl are calculated us- 
ing the empirical fact that correlations between different 
symmetry channels were found to be negligible. Note 
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TABLE VIII. Glueball mass ratios. 



m(2 + 


f )/m(0 ++ ) 


1.39(4) 


m((T 


f )/m(0 ++ ) 


1.50(4) 


m(0 H 


-+)/m(0 ++ ) 


1.54(11) 


m(l+ 


-)/m(0 ++ ) 


1.70(5) 


m(2" 


f )/m(0++) 


1.79(5) 


m(3+ 


-)/m(0++) 


2.06(6) 


m(0*~ 


- + )/m(0 ++ ) 


2.11(6) 


m(0" 


f ) M(2++) 


1.081(12) 



that the pseudoscalar glueball is clearly resolved (at the 
7(7 level) to be heavier than the tensor. 

All of the glueball states shown in Fig. |^ are stable 
against decay to lighter glueballs. In the PC = ++ sec- 
tor, the threshold for decay into two identical ++ glue- 
balls having zero total momentum is at twice the mass of 
the scalar glueball. Although this lies below the mass of 
the 3 ++ glueball, Bose symmetrization prohibits odd L 
partial waves, where L is the relative orbital angular mo- 
mentum, so that the 3 ++ glueball cannot decay into two 
identical scalar glueballs. In the PC = — h sector, the 
lowest-lying two-glueball state consists of the ++ and 
2 ++ glueballs in a relative P-wave; all of our glueballs in 
this sector have masses below the sum of the scalar and 
tensor glueball masses. States of total zero momentum 
and comprised of the ++ and 1^ glueballs with rel- 
ative orbital angular momentum L are the lowest-lying 
two-glueball states in the PC — H — sector when L is 

even and in the PC = sector when L is odd. Only 

the H glueball has sufficient mass to decay into two 
such glueballs; however, this decay is forbidden because 
L — 1 is required to make a state of zero total angular 
momentum. 

To convert our glueball masses into physical units, the 
value of the hadronic scale rg must be specified. We used 
r^ 1 = 410(20) MeV from Ref. §[ to obtain the scale 
shown on the right-hand vertical axis of Fig. ||. This es- 
timate was obtained by combining Wilson action calcu- 
lations of a/ro with values of the lattice spacing a deter- 
mined using quenched simulation results of various phys- 
ical quantities, such as the masses of the p and 4> mesons, 
the decay constant /„•, and the IP— IS* splittings in char- 
monium and bottomonium. Note that the errors shown 
in Fig. H do not include the uncertainty in r^ 1 . For the 
lowest-lying glueballs, we obtain m(0 ++ ) = 1730(50)(80) 
MeV and m(2++) = 2400(25)(120) MeV, where the first 
error comes from the uncertainty in r^mc and the second 
error comes from the uncertainty in 1 . A great deal of 
care should be taken in making direct comparisons with 
experiment since these values neglect the effects of light 
quarks and mixings with nearby conventional mesons. It 
is this mixing which has made the search for an incontro- 
vertible experimental signal so difficult. A glueball hav- 




FIG. 8. The mass spectrum of glueballs in the pure SU(3) 
gauge theory. The masses are given in terms of the hadronic 
scale ro along the left vertical axis and in terms of GeV along 
the right vertical axis (assuming r^ 1 — 410 MeV). The mass 
uncertainties indicated by the vertical extents of the boxes 
do not include the uncertainty in setting ro. The locations 
of states whose interpretation requires further study are indi- 
cated by the dashed hollow boxes. 

ing exotic J will not mix with conventional hadrons 
and would be ideal for establishing the existence of glue- 
balls. Unfortunately, our results indicate that the lightest 
such state, the 2" 1 glueball, has a mass greater than 4 
GeV. 

Kuti has recently pointed out |l^] that the glueball 
spectrum shown in Fig. ^ can be qualitatively understood 
in terms of the interpolating operators of minimal dimen- 
sion which can create glueball states. With the expecta- 
tion that higher dimensional operators create higher mass 
states, the authors in Ref. [^0|, following an approach 
suggested in Refs. constructed all operators of 

dimension 4, 5, and 6 capable of creating glueballs from 
the QCD vacuum. Such operators are gauge-invariant 
combinations of the chromoelectric and chromomagnetic 
fields; operators equivalent to a total derivative or re- 
lated to a conserved current are excluded. The lowest 
dimensional operators capable of creating glueballs are 
of dimension four and have the form Tri^F 1 ^, where 
F^ u is the gauge field strength tensor; these operators 
create glueballs with J PC = 0++,2++,0-+ and 2~+. 



The dimension five operators of the form TrF^DgFap, 
where is the covariant derivative, produce only two 
new glueball states having J PC = 1 ++ and 3 ++ . At 
dimension six, operators of the form T^F^vFapFsa pro- 
duce J PC = ±+ ,l ±± ,2 ±± , and 3 ±_ ; operators of the 
form TTF tlu {D a ,D p }F Sa produce J PC = l-+,3~ + , and 
4 ±+ . Of course, this ordering should not be taken too 
quantitatively, but we find that the method provides a 
reasonably satisfactory explanation of the observed spec- 
trum, especially given the simplicity of the approach. Of 
the lightest six states we resolve, four have the quantum 
numbers expected from the dimension-four interpolating 
operators. The method also explains the absence of any 
low-lying ±_ and 1 h glueballs. 

The spectrum of Fig. || can also be reasonably well 
explained in terms of a simple constituent gluon model 
in which the fundamental gluon field is replaced by the 
Hartree modes of a constituent field with residual per- 
turbative interactions; the Hartree modes are taken to 
be the modes of a free gluon inside a spherical cavity 
with confining boundary conditions. Such a (bag) model 
has been recently revisited in Ref. (ljj. Using values 
for a s and the bag pressure appropriate for heavy-quark 
spectroscopy and the static-quark potential, remarkable 
agreement with the observed levels of Fig. @ was found. 



VII. CONCLUSION 

In this paper, we used numerical simulations of glu- 
ons on spatially-coarse, temporally-fine lattices to signif- 
icantly improve our knowledge of the glueball spectrum 
in SU(3) Yang-Mills theory. This is an important step 
towards understanding glueballs in the real world. Six 
simulations for spatial grid separations ranging from 0.1 
to 0.4 fm were performed on DEC Alpha and Sun Ultra- 
spare workstations. Care was taken to differentiate single 
glueball states from unwanted two-glueball and torelon- 
pair states. An additional small-volume simulation as- 
sisted in the identification of the single glueball states 
and demonstrated the smallness of systematic errors from 
finite volume. The simulation results were extrapolated 
to the continuum limit and the continuum spin quantum 
numbers were identified. The end result, shown in Fig. ||, 
was a nearly complete survey of the glueball spectrum in 
the pure glue theory below 4 GeV. A total of thirteen 
glueballs were found, and two other tentative candidates 
were also located. 

In the future, we plan to improve our determinations 
of the scalar-channel glueballs by simulating with an ac- 
tion that includes an additional two-plaquette interac- 
tion. We also intend to extend our anisotropic lattice 
technology to include quarks. With the help of femto- 
universe techniques, we hope ultimately to investigate 
the properties of glueballs in reality. 
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